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Only some special open surfaces satisfying the shape equation of lipid membranes can be compat- 
ible with the boundary conditions. As a result of this compatibility, the first integral of the shape 
equation should vanish for axisymmetric lipid membranes, from which two theorems of non-existence 
are verified: (i) There is no axisymmetric open membrane being a part of torus satisfying the shape 
equation; (ii) There is no axisymmetric open membrane being a part of a biconcave discodal sur- 
face satisfying the shape equation. Additionally, the shape equation is reduced to a second-order 
differential equation while the boundary conditions are reduced to two equations due to this com- 
patibility. Numerical solutions to the reduced shape equation and boundary conditions agree well 
with the experimental data [A. Saitoh et al, Proc. Natl. Acad. Sci. USA 95, 1026 (1998)]. 

PACS numbers: 87.16.dm, 87.10.Ed 



I. INTRODUCTION 

The elasticity and configuration of lipid vesicles have 
attracted much theoretical attention of physicists [H-Q 
since Helfrich proposed the spontaneous curvature model 
of lipid bilayers in his seminal work Q . The shape equa- 
tion to describe equilibrium configurations of lipid vesi- 
cles was derived in 1987 based on Helfrich's model. 
There are two typical analytical solutions to the shape 
equations. One is a torus with a ratio (-\/2) of its two 
generation radii [1, S. Another is a vesicle with bicon- 
cave discoidal shape [ifl]. In fact, the latter solution does 
not correspond to a vesicle free of external force because 
of a logarithmic singularity in the solution [lTI - [l3l |. 

The opening-up process of lipid vesicles by talin was 
observed by Saitoh et al. 14], which pushes us to in- 
vestigate the shape equation and boundary conditions 
of lipid membranes with free exposed edges. This topic 
was discussed theoretically and numerically by several 
researchers (T5l - [2l| . Based on Helfrich's model, the 
shape equation and boundary conditions were derived by 
Capovilla et al. [H, [ll], Tu et al. [Ulil, and Yin et 
al. [l^ in different forms. Wang and Du obtained var- 
ious shapes of open membranes through numerical sim- 
ulations by phase field method ^2^. Using the area dif- 
ference elasticity model, Umeda et al. derived the shape 
equation and boundary conditions and then compared 
their numerical results with the experiment j2]| . They 
found that the line tension of the free edge of the open 
lipid membrane increases with decreasing the concentra- 
tion of talin |2lj . The above theoretical and numerical re- 
sults can be gene ralized to investigate adhesions between 
lipid vesicles [22, 1231 , configurations of lipid vesicles with 
different lipid domains [3 [13, [H, IHl , and vesicle for- 
mation [2g|. However, the above theoretical researches 
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|15l - ll9j do not contain sufficient discussions on analyti- 
cal solutions to the shape equation with the boundary 
conditions. Additionally, the authors merely compared 
their numerical results with the experimental ones quali- 
tatively in their numerical work i21,] . It is still lack of 
the quantitative comparison between the numerical and 
experimental results. Two natural questions are led to: 
Can we find analytical solutions? At least, it is instruc- 
tive to investigate the possibility of finding the analytical 
solutions. Can we use the numerical results to fit the ex- 
perimental data quantitatively? We hope we can do that 
by taking the number of parameters as small as possible. 

Generally speaking, the shape equation derived from 
Helfrich's model is a fourth-order nonlinear differential 
equation, while the boundary conditions include three 
nonlinear equations describing the shapes of the free 
edges of lipid membranes. In principle, one can obtain 
the general solution with unknown constants to a linear 
differential equation, and then determine the unknown 
constants by using the linear boundary conditions. Thus 
there is no mathematical difficulty to find the solution 
satisfying both the linear differential equation and lin- 
ear boundary conditions. However, the problem becomes 
more complicated if both the differential equation and 
boundary conditions are nonlinear. There is no general 
solution to a nonlinear differential equation in mathe- 
matics. Consequently, one can only conjecture some spe- 
cial solutions in a few cases. If we further consider the 
boundary conditions, only a few ones among the above 
known solutions can fit them. Therefore, it is quite help- 
ful to investigate the feature of the special solutions that 
can satisfy both the nonlinear differential equation and 
the boundary conditions. Since it is very difficult to ob- 
tain solutions to the shape equation with the boundary 
conditions, we may first conjecture a surface satisfying 
the shape equation, and then find a curve in the surface 
satisfying the boundary conditions as an edge of the sur- 
face. However, one might not find any curve satisfying 
the boundary conditions for a given surface satisfying the 



shape equation. Only some special ones among the sur- 
faces satisfying the shape equation can admit the bound- 
ary conditions. The profound reason is that the points in 
the boundary curve should satisfy not only the boundary 
conditions, but also the shape equation because they also 
locate in the surface. In other words, there exist some ad- 
ditional constraints between the shape equation and the 
boundary conditions. These constraints which have not 
been touched in Refs. p"5l - l2]| are called the compatibility 
condition in this paper. 

It is not a straightforward task to find the compati- 
bility condition in general case. The axisymmetric lipid 
membranes with edges will give us some clues. The shape 
equation is reduced to a third-order differential equation 
in axisymmetric case [13, . Zheng and Liu proved that 
it was integrable [1^, and could be further transformed 
into a second-order differential equation with an integral 
constant. In this paper, we will show that the compati- 
bility condition is that this integral constant vanishes for 
axisymmetric membranes. Due to this compatibility, the 
shape equation is reduced to a second-order differential 
equation while the boundary conditions are reduced to 
two equations. The rest of this paper is organized as 
follows: In Sec. |lll we present the general shape equa- 
tion and boundary conditions of lipid membranes with 
free edges. In Sec. IIII[ we discuss the compatibility be- 
tween the shape equation and boundary conditions in 
axisymmetric case, and then verify two theorems of non- 
existence. In Sec lIVl we find some axisymmetric numer- 
ical solutions and compare them with experimental data 
quantitatively. A brief summary is given in the last sec- 
tion. 




FIG. 1: Right-handed frame of an open surface with a bound- 
ary curve C. 



spontaneous curvature. A and 7 are the surface tension 
and line tension, respectively. 

By calculating the variation of free energy ([T]) , we can 
obtain [l3| 

{2H + co)(2if2 - cqH - 2K) - 2~XH + V^(2H) = 0, (2) 



and 



C 



[{2H + Co 
[~2dH /de-z + ■ynn + kdrg/ds] 



0, 



[{1/2){2H + co)2 + kK+~X + ^Kg 



= 0, 



(3) 
(4) 
(5) 



II. SHAPE EQUATION AND BOUNDARY 
CONDITIONS 



A lipid membrane with a free edge is represented as 
an open surface with a boundary curve C. As shown 
in Fig. [U we can construct an orthogonal right-handed 
frame {ei, 62, 63} at each point of the surface such that 
83 is the normal vector of the surface. For each point in 
the boundary curve, we take 62 to be perpendicular to 
the tangent direction ei of the boundary curve and point 
at the side that the surface is located. 

The free energy of the membrane may be expressed as 



F 



[{kc/2){2H + cof + kK]dA + XA + -fL, (1) 



where the first and second terms are the bending energy 
and the surface energy of the membrane, respectively, 
while the third term is the line energy of the free exposed 
edge. H and K are the mean curvature and gaussian 
curvature of the surface, respectively. dA is the area 
element of the surface. A and L are the total area of 
the surface and the total length of the boundary curve, 
respectively, kc and k are the bending moduli, cq is the 



where A = A/fcc, k = k/kc, 7 = j/kc are the reduced 
surface tension, reduced bending modulus, and reduced 
line tension, respectively. k„, Kg, t^, and ds are the 
normal curvature, geodesic curvature, geodesic torsion, 
and arc length element of the boundary curve, respec- 
tively. Equation ([2]) determines the equilibrium shape of 
the membrane, thus we call it shape equation. For a given 
surface satisfying the shape equation, Eqs. 131)-® deter- 
mine the shape of the boundary curve and its position 
in the surface, thus we call them boundary conditions. 
Equation ^ expresses the normal force balance of the 
membrane. Equation ^ is the moment balance equa- 
tion around ei at each point in curve C. Equations ^ 
and (O are the force balance equations along 63 and 62 
at each point in curve C, respectively. Thus, in general, 
the above four equations are independent of each other. 



III. COMPATIBILITY BETWEEN THE SHAPE 
EQUATION AND BOUNDARY CONDITIONS 

We have mentioned that only some special ones among 
the surfaces satisfying the shape equation ^ can admit 
the boundary conditions ([S])-®. What is the common 
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FIG. 2: Outline of an open surface. Each open surface can be 
generated by a planar curve AC rotating around z axis, i/' is 
the angle between the tangent line and the horizontal plane. 



boundary conditions (O-®. In particular, the points in 
the boundary curve should satisfy not only the boundary 
conditions, but also the shape equation (fTO|) because they 
also locate in the surface. That is, Eqs. ©-(HI) and ([TU]) 
should be compatible with each other in the edge. Substi- 
tuting Eqs. (Ill)-® into (fTO|). we derive the compatibility 
condition to be 



Vo=0- 



(11) 



No we will discuss two examples and verify two theo- 
rems of non-existence by considering the above compat- 
ibility condition. 




feature of these special surfaces? we will find this feature 
for axisymmetric surfaces. 

When a planar curve AC shown in Figl2] revolves 
around z axis, an axisymmetric surface is generated. Let 
^p represent the angle between the tangent line and the 
horizontal plane. Each point in the surface can be ex- 
pressed as vector form r = {pcos^, p sin cj), z(p)} where p 
and (j) are radius and azimuth angle that the point corre- 
sponds to. Introduce a notation a such that cr = 1 if ei 
is parallel to dr/dcj), and cr = —1 if ei is antiparallel to 
dvjdc^) in the boundary curve generated by point C. The 
above equations Q-® are transformed into 



- Co) — + 



2 



cph 
2 



2 



2K 



Xh 



COS'0 



{pcoail^h'Y = 0,(6) 



h-co + fcsin ijj/p =0,(7) 



-acos^ph' -f 7sin V'/pIc — 0,(8) 



simp 



kK + X 



_ COS tp 
CT7 

p 



- 0,(9) 



c 



where h = siTi^p/p+ (sin^)' and K = sin-ip^simp)' / p. 
The 'prime' represents the derivative with respect to p. 

The shape equation ([6]) is a third-order differential 
equation. Following Zheng and Liu's work j29j . we can 
transform it into a second order differential equation 



costph' 



{h — Co) simpip' 
rjo tan ip 
p cos ip 2 



— X tan ^p 



(10) 



with an integral constant r/o (so called the first integral). 
The configuration of an axisymmetric open lipid mem- 
brane should satisfy the shape equation ([5]) or ([TU)) and 



FIG. 3: (Color online) Two non-existent axisymmetric open 
membranes: (a) A part of a torus; (b) A part of biconcave 
discodal surface. 

First, let us consider a part of a torus shown in Fig. [3^ 
generated by an arc expressed by sin ^ = ap+p^ with two 
non- vanishing constants a and /3. Substituting it into the 
shape equation (fTOl) . we obtain co = 0, /3 = ^/2, A = 0, 
and ?7o = ~a. That is, the torus can be a solution to the 
shape equation. However, rjQ = —a ^ contradicts to 
the compatibility condition (|lip . Thus we arrive at: 

Theorem 1. There is no axisymmetric open membrane 
being a part of torus generated by a circle expressed by 
sin-i/) = ap + \f2. 

Secondly, we consider a biconcave discodal surface 
[loj generated by a planar curve expressed by smip = 
ap\'n.{p/P3) with two non- vanishing constants a and /3. 
To avoid the logarithmic singularity at two poles, we may 
dig two holes around the poles in the surface as shown in 
Fig. Eb- Substituting sinip = ap\n{p/pi) into the shape 
equation (fTO|) . we obtain A = 0, a = cq, and ryo = — 2co. 
That is, the biconcave discodal surface can be a solution 
to the shape equation. However, rjQ = —2cq ^ contra- 
dicts to the compatibility condition (fTTj) . Thus we arrive 
at: 

Theorem 2. There is no axisymmetric open membrane 
being a part of a biconcave discodal surface generated by 
a planar curve expressed by simp = co/9ln(p//3). 

In the above discussion, the theorems of non-existence 
are deduced as natural corollaries of the compatibility 
condition. It does not mean that the proofs are unique. 
The other proofs are presented in Appendix [X] In 
Ref. [131, the present author has proved that there is 
no open lipid membrane being a part of a constant mean 
curvature surface. These theorems reveal that it is al- 
most hopeless to find analytical solutions to the shape 
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equations with the boundary conditions. Thus we need 
to seek for numerical solutions. 



IV. AXISYMMETRIC NUMERICAL 
SOLUTIONS 

The compatibility condition leads to a more important 
result that the shape equation can be simplified as 



3.0 ■ 



2.5 ■ 



cos iph' + (/i — Co) sin ip^p' — X tan ip - 



tani/> 



while the boundary conditions can be taken only two 
equations ([7]) and ^ because Eq. ([SJ is not independent 
of Eqs. (IZl), (12), and However, it is still very difficult 
to obtain analytical solutions to Eq. (jl2p with boundary 
conditions ([7]) and ([9]). We will find axisymmetric nu- 
merical solutions and compare them with experimental 
data [1^ in this section. 

Because "0 might be a multi-valued function of the in- 
dependent variable p, the above equations are unsuitable 
for numerical solutions. Here we take the arc-length of 
curve AC in Fig. [2] as an independent variable. Then we 
have p = cosip and z = sin?/;, where the 'dot' represents 
the derivative with respect to the arc-length. The shape 
equation can be transformed into 



^ = - 



tan -0 • 2 cos ipip sin 2'ip 



-r- 



p 



2p2 



-Atani/)- 



while the boundary conditions become 

sinVj 



tani/) f sin0 
P 

(13) 



V' - Co + 



(l-Hfc) 



P 



0, 



and 



fcco 



sin0 
P 




sin?/' 



+ A + 7 



cos-0 

p 



(14) 



= 0. 



(15) 

In fact, these equations can be also derived from the La- 
grangian method as shown in Appendix |B] In addition, 
we impose the initial conditions 2(0) — p{Q) — /xm, and 
?/'(0) — 0. We can use the shooting method to find nu- 
merical solutions to Eq. ()13p with boundary conditions 
[Eqs. (fn|) and (IT5|) ] and these initial conditions, and then 
fit the parameters (fc, Co, A, 7) with experimental data. 
The basic idea is as below. For the given values of fc, 
Co, A, 7 and ■0(0), we can solve Eq. (fT3|) with boundary 
conditions (|14l) and (flSI) . Then we compare the graph of 
the solution to the outline of the open membrane in the 
experiment. Tune the values of the parameters until the 
graph of the solution and the outline of the membrane 
almost superpose each other. Thus we obtain a group of 
proper values of the par ameters (fc, cq. A, 7). 

In the experiment |14| , the hole of the lipid membrane 
is enlarged with increasing the concentration of talin, and 
vice versa. Talin molecules adhere to the edge of the 
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FIG. 4: (Color online) Numerical results (solid, dash, and dot 
lines) and experimental data (squares, circles, and triangles 
extracted from Fig. 3.1 to K in Ref. [ij]) of the outlines of 
an axisymmetric lipid membrane at different concentration of 
talin. 



membrane. Thus it is reasonable to assume that the line 
tension of the edge depends on the concentration of talin, 
while the bending moduli and spontaneous curvature of 
me membrane do not. That is, we should have the com- 
moii values of k and Cq for a membrane at different con- 
centrations of talin. This gives a constraint in our fitting. 
As shown in Fig. 21 our numerical results (solid, dash, and 
dot lines) obtained from Eqs. (IT3)) - (IT5t agree well with 
the experimental data (squares, circles, and triangles ex- 
tracted from the outlines of the membrane with decreas- 
ing the concentration of talin) . The common parameters 
are fitted as fc = —0.122 and co = 0.4 /im^^. The neg- 
ative k reveals that the surface with positive Gaussian 
curvature is more energetically favorable than that with 
negative Gaussian curvature. The positive co reflects the 
asymmetry of the bilayer composed in the experiment 
[14l|, which makes the membrane bend like a standing 
upright cup. The other parameters are shown in the fig- 
ure. The reduced line tension 7 increases from 0.66 pm~^ 
to 0.79 /zm~^ with decreasing the concentration of talin, 
which is the same as the conclusion of Ref. [2l|. The 
surface tension depends on the shapes and line tension. 
Intuitively, the line tension of the edge induces compres- 
sion stress in the membranes with the similar shapes gen- 
erated by solid line in Fig. |4l thus the surface tension is 
negative. On the contrary, the tension of the edge in- 
duces stretching stress in the membranes with the similar 
shapes generated by dash or dot lines in Fig. 21 thus the 
surface tension is positive. The variation of surface ten- 
sion can be comprehensively understood from Eq. (jA4p . 
For the membrane in Fig.|4l 2H = —{sinip/p+ip) < and 
decreases according to the sequence of the solid, dash and 
dot lines. Therefore, the surface tension of the membrane 
increases according to the same sequence with consider- 
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ing Eq. (jA4[) . Furthermore, we also examine that our 
numerical results indeed satisfy the constraint (|A4p . 



V. CONCLUSION 

In the above discussion, we investigate the compati- 
bility between shape equation and boundary conditions 
of lipid membranes with free edges. The main results 
obtained in this paper are as follows. 

(i) The compatibility condition for axisymmetric lipid 
membranes with free edges is that the first integral of the 
shape equation © should be vanishing, i.e., Eq. ([TT|). 

(ii) Two theorems ( Theorem 1 and Theorem 2 in 
Sec. IIII[) of non-existence are verified as natural corol- 
laries of the compatibility condition, which give two ex- 
amples to reveal that one indeed might not find any curve 
satisfying the boundary conditions in a given surface sat- 
isfying the shape equation. These theorems also correct 
two flaws on analytical solutions in Ref. [l7j . 

(iii) The shape equation of axisymmetric lipid mem- 
branes is reduced to Eq. (fT2|) . Then only two equations 
in boundary conditions are independent. This conclusion 
is the same as the case in Ref. [l^ with vanishing k and 
Co- 

(iv) As shown in Fig. 21 the numerical solutions to the 
reduced shape equation (IT2]) with boundary conditions 
([7]) and (0) agree well with the experimental data . 

Finally, we would like to point out two difficulties that 
we have not fully overcome yet: (i) The compatibility 
condition between shape equation and boundary con- 
ditions for asymmetric (not axisymmetric) lipid mem- 
branes with edges is unclear. We do not even know 
whether it exists, much less what it is. (ii) We use 
the shooting method to find numerical solutions. But 
this method is not so efficient to the numerical solutions 
due to the complicated boundary conditions. A much 
more efficient method is expected. The above challenges 
should be addressed in the future work. 
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Appendix A: Other proofs to theorems 
non-existence 

The proofs can be divided into two classes in terms 
of different starting points. One is based on the stress 
analysis, another is based on the scaling argument. 



1. Stress analysis 

Capovilla et al. proposed the stress tensor in a lipid 
membrane and then derived the shape equation and 
boundary conditions from the stress tensor [1^, [l^ . Re- 
cently, they found that [H, Ell the line integral 

j> dsl^ia -z^c, (Al) 

where F is any circle perpendicular to the symmetric axis 
in an axisymmetric membrane. 1°" and fa represent the 
normal of F tangent to the membrane surface and the 
stress in the membrane, respectively, z is the unit vector 
along the symmetric axis, c is a constant dependent on 
the topology and the curvature singularity of the mem- 
brane. 

First, the constant c is non-vanishing for an axisym- 
metric torus free of curvature singularity, which implies 
that the stress in each circle F perpendicular to the sym- 
metric axis in the torus surface cannot be zero. However, 
the stress in the free edges should be vanishing. Thus, 
we cannot find any F as a free edge of an axisymmetric 
open membrane being a part of the torus, i.e., theorem 
1 is arrived at. 

Secondly, there exist singularity points at two poles 
of the biconcave discodal surface generated by a planar 
curve expressed by sin?/' = oip\p.(pl 0). The singularity 
results in a non- vanishing c [12l . [l3| , and then non-zero 
stress in each circle F in the biconcave discodal surface. 
Thus, we cannot find any F as a free edge of an axisym- 
metric open membrane being a part of the biconcave dis- 
codal surface, i.e., theorem 2 is arrived at. 

2. Scaling argument 

The free energy ([T]) can be written in another form 
F = j[{kc/2){2Hf + kK]dA 

+ 2kcCo J HdA + {X + kccl/2)A + -/L. (A2) 

Let us consider the scaling transformation r — !■ Ar, where 
the vector r represents the position of each point in 
the membrane and A is a scaling parameter |16|. Un- 
der this transformation, we have A -> A^A, L — > AL, 
H A-^H, and K A-^K. Thus, Eq. ^M^, is trans- 
formed into 

F(A) = J[{kc/2){2Hf + kK]dA 

+ 2kcCoA J HdA+{X + kccl/2)A^A + jA(JA3) 

The equilibrium configuration should satisfy dF/dA — 
when A = 1 [Tg]. Thus we obtain 

2co / HdA + (2A + cl)A + jL^O. (A4) 
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This equation is an additional constraint for open mem- 
branes. 

As shown in Sec. Illli if there exists an open membrane 
being a part of torus, then the shape equation ^ requires 
A = and co = 0, which contradicts the constraint (|A4p 
because 7L > 0. Thus we arrive at theorem 1. 

Because WiUmore surfaces satisfy the special form of 
Eq. ^ with vanishing A and Cq [31[, as a byproduct of 
the constraint (jA4l) . we obtain a much stronger theorem 
of non-existence: There is no open membrane being a 
part of a Willmore surface. 

Next, we turn to the biconcave discodal surface. If 
there exists an open membrane being a part of a bicon- 
cave discodal surface generated by a planar curve ex- 
pressed by sin-!/; = copln(p//3), the shape equation ([2]) 
requires A = 0. Substituting 2H = -co[l + 21n(p//3)] 
into Eq. (|A4p . we will not obtain a contradiction. Thus 
theorem 2 cannot be deduced from the scaling argument. 



Jo \dp ds dp) ^ ^ c»p ^ 

n / d dC\ ^ ^ dC^ 
+ / 0^1-1^ '^^c?s+ —Sz 
Jq \ ds dz J 

+ / {p — COS -tp) d(^ds + / {z — sin ^) Sr]ds 



Jo Jo 

- n\cSs2 = 0, (B3) 

where the Hamiltonian = + + - ^- Im- 
posing ip{0) — 0, p(0) — z(0) — /im, and substituting 
Eq. ([§21 into Eq. (|Mt . we can obtain 

CsinV-p/ = 0, (B4) 

77 — constant, (B5) 

p — costp, (B6) 

i; = sinV', (B7) 



Appendix B: Derivation of the reduced shape 
equation and boundary conditions by using the 
Lagrange method 

For the revolving surface generated by the planar curve 
shown in Fig. [21 Eq. ([ij can be transformed into 

F/2TTkc= [pf/2 + ksmil;ij + Xp + ^p\ds, (Bl) 
Jo 

with / = sin?/'//? + ip ~ cq. We should minimize F/2'Kkc 
with the constraints p — cosip and i = sin'0, thus we 
construct an action S — J^^ Cds with a Lagrangian ^ 

C = pf'^ /2 + ksmip^ + Xp+jp+C{p—costp)+r]{z — sin^), 

(B2) 

where C and rj are two Lagrange multipliers. In terms of 
the variational theory, we can derive 



SS = 



with boundary conditions 



SCds — 'Hds2 
/o \diP dsdipj dip 



f 



~ s'mtb 
k 



0, 



Clc + 7 = 0, 

n\c = o. 



(B8) 

(B9) 
(BIO) 
(Bll) 



Eq. (jBSp is equivalent to boundary condition 

Because C does not explicitly contain s, 'H is a con- 
stant. Combining Eqs. (|55|) . (|B10|) . (|Blip and the defi- 
nition of we derive 



CcosV-p[/(/-27A)/2 + A] 



(B12) 



From Eqs. (|B4p and (|B12I) we can obtain the shape equa- 
tion (113]). Equation (IBTT|) can be transformed into the 
boundary condition p3|) with Eq. (jB8|) . From Eqs. (jB4[) 
and (|B9[) we can also obtain the other boundary condi- 
tion, which is not independent of the shape equation (jl3p 
and boundary conditions (ITil) and 
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